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The theory of the magnetic catalysis of chiral symmetry breaking in QED is developed. An 
approximation for the Schwinger-Dyson equations describing reliably this phenomenon is established, 
i.e., it is shown that there exists a consistent truncation of those equations in this problem. The 
equations are solved both analytically and numerically, and the dynamical mass of fermions is 
determined. 
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I. INTRODUCTION 

Recently the magnetic catalysis of dynamical chiral symmetry breaking has been established as a universal phe- 
nomenon in 2+1 and 3+1 dimensions: a constant magnetic field leads to the generation of a fermion dynamical 
mass even at the weakest attractive interaction between fermions ||-||] . The essence of this effect is the dimensional 
reduction D ^ D — 2 in the dynamics of fermion pairing in a magnetic field: at weak coupling, this dynamics is 
dominated by the lowest Landau level (LLL) which is essentially {D — 2)-dimensional ||l|-^. The effect may have 
interesting applications in condensed matter physics and cosmology |0,||-B|. 

In particular, this phenomenon was considered in 3+1 dimensional QED [7|,p|~p^. Since the dynamics of the 

LLL is long-range (infrared), and the QED coupling constant is weak in the infrared region, one may think that the 
rainbow (ladder) approximation is reliable in this problem. As was shown in Refs. (|,^,^, in all the covariant gauges, 
the dynamical mass of fermions in this approximation is 



CVI^exp 



TT / TT \ 

2 V2^/ 



1/2- 



(1) 



where B is a magnetic field, the constant C is of order one and a is the renormalized coupling constant related to the 
scale /x^ '-^ |ei?|. 

Are higher order contributions indeed suppressed in this problem? The answer is "no". As was shown in Refs. [^D, 
because of the (l+l)-dimensional form of the fermion propagator of the LLL fermions, there are relevant higher order 
contributions. In particular, considering this problem in the improved rainbow approximation (with the bare vertex 
in the Schwinger-Dyson equations for both the fermion propagator and the polarization operator), it was shown that, 
in all the covariant gauges, the fermion mass mdyn is given by Eq. (|l|) but with a ^ a/2 

As we wrote in the paper ^ , "it is a challenge to define the class of all those diagrams in QED in a magnetic field 
that give a relevant contribution in this problem". In this paper, we will solve the problem. (A brief outline of our 
results was given in Ref. [p^.) We will show that there exists a (non-covariant) gauge in which the Schwinger-Dyson 
equations written in the improved rainbow approximation are reliable: in other words, in that gauge, there exists a 
consistent truncation of the Schwinger-Dyson equations for this TioTi-pGTtuTbdtivc problem. The expression for uifj^yfi 
takes the following form, 



exp 



a\n{Ci/Na) 



(2) 



where N is the number of fermion flavors, i^(a) ~ {NaY^^ , Ci ~ 1.82 ± 0.06 and the constant C is of order one. 

This expression for mdyn is essentially different from that in the rainbow approximation (Q). As we will see, this 
reflects rather rich and sophisticated dynamics in this problem. 



*On leave of absence from Bogolyubov Institute for Theoretical Physics, 252143, Kiev, Ukraine. 
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The paper is organized as follows. In Section || the Schwinger-Dyson equations in QED in a magnetic field are 
discussed. In Section III we define a (non-covariant) gauge in which the improved rainbow approximation for these 
equations is reliable. In Section IV the loop expansion for the Schwinger-Dyson equations is considered and it is 
confirmed that the improved rainbow approximation is indeed reliable in the gauge introduced in Section 



III, 



In 

Section M the Schwinger-Dyson equations in this approximation are solved both analytically and numerically. In 
Section |V]| we summarize the main results of the paper. The details of our analysis and some useful formulas and 



relations are presented in Appendices ^ and 



II. THE SCHWINGER-DYSON EQUATIONS IN QED IN A MAGNETIC FIELD. 

The Lagrangian density of massless QED in a magnetic field is 

>C = --^F^'^F^, + i {ii^D,.)i,] + J^A^, (3) 



where the covariant derivative is 



D^ = d^- ieA^ (4) 



and the source provides a chosen external field: (0|Ap|0) = A'^f*. As we will see below, in the case of a constant 
magnetic field B, with 

(0'-f^2,fxi,o), (5) 

the source J'' = 0. Notice that the vector potential A^^* (^) corresponds to the so called symmetric gauge for the 
external vector potential, and the magnetic field is in the -1-2:3 direction. 

Besides the Dirac index (n), the fermion field carries an additional flavor index a = 1,2, . . . , N. Then the Lagrangian 
density in Eq. (0) is invariant under the chiral SUl{N) x SUji{N) x Uv{l) symmetry (we will discuss the anomalous 
Ua{1) in Sec. |V|). 

The Schwinger-Dyson (SD) equations in QED in external fields were derived by Schwinger and Fradkin (for a 
review, see Ref. The equations for the fermion propagator G{x,y) are 

G-\x,y) = S-\x,y) + J:{x,y), (6) 

T.{x, y) = 47ra7'' / G{x, z)T''{z, y, z')V^f,{z', x)d*zd*z'. (7) 



Here S{x,y) is the bare fermion propagator in the external field A^^*, E(x,y) is the fermion mass operator, and 
T^ij.u{x, y), r'^(x, y, z) are the full photon propagator and the full amputated vertex. 
The full photon propagator satisfies the equations 

V;^ix, y) D-l{x -y) + Il,,,{x, y), (8) 

n^i/(a;, y) = -47rQ:tr7p j d'^ud'^zG{x,u)Ty{u,z,y)G{z,x), (9) 



where D^^{x — y) is the free photon propagator and 11^1,(2;, y) is the polarization operator. 
The equation for the external photon field A^^^ = (0|A^|0) in a covariant gauge is 

uAlf - Xd^d'-AT = -J^ - (OIjJO), (10) 

where the vacuum current (0|jp|0) is (0|j^|0) — — etr(7^G'(a;, x)) and A is the gauge parameter. Notice that for A^^* 
in Eq. (^), corresponding to a constant magnetic field, the source is zero. Indeed, in this case, = — (0|j;i|0), and 
(0|j^|0) = -etr(7^G(a;, x)) = because of the symmetry S0{2) x SO{l, 1), with 50(2) and SO{l, 1) corresponding 
to rotations in the xi — a;2-plane and Lorentz transformations in the xq — xa-hyperplane, respectively. 

The bare fermion propagator S{x, y) in a constant magnetic field was calculated by Schwinger [l^ . In the symmetric 
gauge (^), it has the form 

S{x, y) = exp {zex^^Al^'iy)) S{x - y), (11) 
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where the Fourier transform of S is 



oo 

S{p) ~ J dsexp 





2 2 2 tanieBs) 



eBs 

■ [(pV -P^7^ + +7^7^ tan(eBs)) -p_L7±(l + tan2(eBs))] . (12) 

Here a transverse vector p± = [p^ ,p'^) and m is a fermion mass. Then, using the identity 

ix^Af[z) + iz^Af{y) = ix^A^-^y) + t{x ~ yTAl^^z - y) (13) 

for the vector potential (|^), it is not difficult to show directly from the SD equations (^), (Q), (||) and (||) that 

G{x, y) = exp {zex^Al^^y)) G{x - y), (14a) 

r(x, y, z) = exp [iex^^Al^'iy)) f{x -z,y-z), (14b) 

'Dt,„{x,y) ='Df,u{x-y), (14c) 

Uf,^{x,y) ^Uf,^{x- y). (14d) 

In other words, in a constant magnetic field, the Schwinger phase is universal for Green functions containing one 
fermion field, one antifcrmion field, and any number of photon fields, and the full photon propagator is translation 
invariant. 

Our aim is to show that there exists a gauge in which the approximation with a bare vertex, 

T''ix,y,z) = j^5{x-y)6ix-z), (15) 
is reliable for the description of spontaneous chiral symmetry breaking in a magnetic field. 

III. NON-COVARIANT GAUGE AND THE IMPROVED RAINBOW APPROXIMATION FOR THE SD 

EQUATIONS 

In this section, we will show that there is a (non-covariant) gauge in which the approximation with a bare vertex 
( |l5| ) (the improved rainbow approximation) is reliable. 

We begin by recalling the following facts concerning the problem of the magnetic catalysis of chiral symmetry 
breaking ||]: 

1. At weak coupling, there is the LLL dominance in the dynamics of fermion pairing. It is because of the presence 
of the large Landau gap of order -y/|e_B|, which is much larger than the dynamical fermion mass mdyn (for weak 
coupling). In other words, higher Landau levels decouple from the infrared dynamics (with k •\/|ei?|) in 
the same way as it happens with Kaluza-Klein (KK) modes in KK theories of gravity. This fact was explicitly 
shown in the Nambu-Jona-Lasinio model and in QED [p^ . 

2. The propagator S{p) ( |l^ ) can be expanded over the Landau levels [p^P]. The contribution from the LLL is 

/ ,\2 Pn + m , ^ 

SLLLip)^2ie~^P^'^ ^ ^O(-), (16) 

Pj| — TO^ 

where the magnetic length I = |ei?|~^/^, p± = (p^,p^), p\\ = {p^,p^), and j5|| = ^"7*' — ^^7"^. The matrix 
O(-) = (1 - n^jhign{eB)) /2 is the projection operator on the fermion states with the spin polarized along 
the magnetic field. This point and Eq. ([iq ) clearly reflect the (l+l)-dimensional character of the dynamics of 
fermions in the LLL. This property is preserved also in the case when the fermion mass is generated dynamically 
(m = mdyn) Q- 



^As was shown in Refs. plpOt|, despite the dimensional reduction 3+1^1 + 1 in the fermion propagator in a magnetic field. 
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3. The presence of the projection operator ^ imphes that the bare vertex for fermions from the LLL is 

= 0^~^7jj'. Therefore the LLL fermions couple only to the longitudinal (0,3) components of 
the photon field. 

4. In the one-loop approximation, with fermions from the LLL, the photon propagator takes the following form in 
covariant gauges P|, 




A 



(17) 



where the symbols _L and || in g^j^i, are related to the (1,2) and (0,3) components, respectively, and A is a 
gauge parameter. The explicit expression for Xl[q\,q^^) = cxp[— ((7^Z)^/2]n(gp is given in Refs. |22|J^. For our 
purposes, it is sufficient to know its asymptotes, 

n(gp-:^^, as |<7i||«m2^„, (18) 

9n 2a \eB\ ,9,9 , 

n(gf)^ V as k^|»m^^„, (19) 

where a = Na. 

Notice that the polarization effects are absent in the transverse components of 'D^i,{q). This is because, as was 
already pointed out in item 3 above, fermions from the LLL couple only to the longitudinal components of the 
photon field. 

5. Then, there is a strong screening effect in the (^g^^u — 9/1 component of the photon propagator. Eq. ( p^ ) 
implies that 

(20) 



q^ + qlU{ql,ql) " q^ - M^' 
with 



= -^\eB\ (21) 



TT 



for ^ |q^| ^ \eB\ and \q'j^\ ^ \eB\. This is reminiscent of the Higgs effect in the (l+l)-dimensional QED 
(Schwinger model) ||,||. 

We emphasize that the infrared dynamics in this problem is very different from that in the Schwinger model: since 
photon is neutral, there is no dimensional reduction for its field in a magnetic field, and there is the four-dimensional 
= — in the photon propagator (pO|). However, the tensor and the spinor structure of this dynamics are exactly 

the same as in the Schwinger model: indeed, the LLL fermion propagator ( p^ and the vertex 0'-~^7^0'-~-' = 0^~'^j^ 

are two-dimensional, and only the longitudinal (0, 3) components of a photon field are relevant here. This point will 
be crucial for finding a gauge in which the improved rainbow approximation (with the bare vertex (p^) is reliable|^. 



there is a genuine spontaneous chiral symmetry breaking in this problem: the Mermin- Wagner- Coleman (MWC) theorem 
forbidding spontaneous breakdown of continuous symmetries at D = 1 -I- 1, is not applicable to this case. The point is that 
the MWC theorem is based on the fact that gapless Nambu-Goldstone (NG) bosons cannot exist in 1 -I- 1 dimensions. On the 
other hand, since the NG bosons connected with spontaneous chiral symmetry breaking are neutral, there is no dimensional 
reduction of their propagator in a magnetic field ]3[|2C|]. Therefore their propagator is (3 + l)-dimensional and there are no 
obstacles for their existence in this case. 

^Since an external magnetic field does not lead to confinement of fermions, their mass is gauge invariant in QED in a magnetic 
field. Therefore any gauge can be used for the calculations of the mass if either the calculations provide the exact result or a 
good approximation is used: i.e., one can show that corrections to the obtained result are small. Below we will define such a 
gauge in this model. 
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We recall that, as was shown in Ref. despite the smallness of a, the expansion in a is broken in covariant gauges 
in this problem. The reason is that, because of the smallness of uidyn in Eq. ([|) as compared to y^eSf, there are 
mass singularities, \ii[\eB\/m?^y^) ^ a"^/^, in infrared dynamics. In Appendix we analyze these singularities in the 
vertex. Calculating the one- loop correction to the vertex, one finds that, when external momenta are of order m^yn 

or less, there are contributions of order aln^ (leBl/m^y^J ~ They come from the term q^^ql / q'^ q'^^ in T>^^{q) in 

Eq. ®. 

How can one avoid such mass singularities? A solution is suggested by the Schwinger model. It is known that there 
is a gauge in which the full vertex is just the bare one It is the gauge with a bare photon propagator 

D^,{k) = -.1 [g^, - ^) - (22) 

with the (non-local) gauge function d = l/(l-fn), where the polarization function n(fc^) = — e^/7rfc^ in the Schwinger 
model (of course, here a, /3 = 0, 1). Then, the full propagator is proportional to gap, 

Va,{k) = Dap{k) + ^ i^9a, - j fc2(i + n(fc^)) ^ -'m^TWm)- ^ ^ 

The point is that since now Vapik) ^ gap and since the fcrmion mass m = in the Schwinger model, all loop 
contributions to the vertex are proportional to 

P2„+i=7a7Ai...7Wi7" = (24) 

in that gauge and, therefore, disappear^. 

Let us return to the present problem. As it was emphasized above, the tensor and the spinor structure of the LLL 
dynamics is (l+l)-dimcnsional. Now, take the bare propagator 



r V T ) " q^q\ 

with d = — g|n/[q^ + g^II] + q^^/q^ ■ Then, the full propagator is 



- DM + * \ „2 I ^2[g2 + 92n(^2 ^ g2)] 




qln{ql,q^^ 



.qjiqv +qj^ql + qlqu 
Wy • ^''^ 

The crucial point is that, as was pointed out above, the transverse degrees of freedom decouple from the LLL dynamics. 
Therefore only the first term in ^^^^{q), proportional to yjlj/, is relevant. 

Notice now that dangerous mass singularities in loop corrections to the vertex might potentially occur only in the 
terms containing gj' = ql^'^ — qf"f^ from a numerator (q| -|- rridyn) of each fermion propagator in a diagram (all other 
terms contain positive powers of mdyn, coming from at least some of the numerators and, therefore, are harmless). 
However, because of the same reasons as in the gauge (^3|) in the Schwinger model, all those potentially dangerous 
terms disappear in the gauge (p^). Therefore all the loop corrections to the vertex are suppressed by positive powers of 
a in this gauge. This in turn implies that those loop corrections may result only in a change C ^ 0(1) C ^ 0(1) 
in expression (|^.) In other words, in gauge (|2^ ) there exists a consistent truncation of the SD equations and the 
problem is essentially soluble in that gauge. ^ 

In the next section, we will consider the loop expansion for the SD equations in this problem in more detail. 



^^"271+1 = follows from the two identities for the two-dimensional Dirac matrices: 7a7A7" ~ and "fx-jx-^-^ — (?AiAi_|_i -I- 
£\i\i+il5 (75 = 7o7i> Sap = -epc £01 = !)• 

''The gauge is unique in that. In other gauges, there is an infinite set of diagrams giving relevant contributions to the 
vertex. Therefore, in other gauges, one needs to sum up an infinite set of diagrams to recover the same result for the fermion 
mass. 
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IV. THE LOOP EXPANSION FOR THE SD EQUATIONS 



The consideration of mass singularities in loop corrections given at the end of the previous section was, though 
general, somewhat heuristic. First of all, one has to define more rigorously the perturbative expansion for the SD 
equations which is used in this problem. It is the loopexpansion based on the CJT (Cornwall, Jackiw, and Tomboulis) 
effective action r(G, P^j/) for composite operators (for a review see Ref. |26|). The conditions for extrema of T 
yield the SD equations, 



6r 



SG{x,y) 



= 0, 



ST 



= 0. 



(27) 



In the loop expansion for F, the full photon and fcrmion propagators are used in two-particle irreducible diagrams for 
F. In QED, the problem is essentially reduced to the loop expansion (with the full photon and fermion propagators) 
for the vertex. 

The full photon propagator is given by Eq. (26), and the full propagator for fcrmions from the LLL has the form 



G{p) = 2ie-(f-L')' 



aHpI)p1^b^{pI) 



(28) 



(compare with Eq. (16) and see below). Here S(p|) is a dynamical mass function of fermions. 

In this section we will derive the SD equations for the fermion propagator in the one-loop and two-loop approxima- 
tions. It will be shown, that while the one-loop approximation, coinciding with the improved rainbow approximation, 
is reliable in non-covariant gauge (p6|), it is not reliable in covariant gauges (p^). 

From Eqs. (^) and (^ one gets the following equation for the fermion prop agat or G{x,y) in the two-loop approxi- 
mation (by using the vertex function in the one- loop approximation, see Eq. (14b) and Eq. ( |Al| ) in Appendix " 

G{x, y) = S{x, y) - Ana / (fxid^yiSix, xi)^''G{xi,yi)^''G{yi,y)'Df,^{xi - yi) 



+ (47ra)2 J d^xid\2d^yid%2S{x, xih>'G{xi,X2h''G{x2,y2h''G{y2,yih''G{yi,y)V^4xi - y2)V,p{x2 - yi). (29) 
Here S{x, y) is the bare fermion propagator of massless fermions (m = 0). The graphic form of this equation is shown 



.0. 



in Fig. pi. After extracting the Schwinger phase factors in the full and bare fermion propagators [see Eq. (14a)], 



G{x,y) 



G{x~y), Six,y)^e^-"<^'^y^Six-y), 



(30) 



Eq. ( p9| ) reads as 

G{x) = S{x) - Ana J d''xid''yie"^^^'^+'^'^^y'^ S{x - xi)rG{x, - yi)rG{yi)Vp,{xi - yi) 

xSix - xi^f^Gixi - X2h^Gix2 - y2h"Giy2 - yi)YG{yi)Vp,{xi - y2)V,p{x2 - yi) (31) 

where A^^* is given in Eq. (|^) and the shorthand xA'^^*{y) stands for a;^Ajj^*(?/). 

First, let us show that the solution to the above equation, G{x), allows the factorization of the dependence on the 
parallel and perpendicular coordinates. 



G{x)^ 



2nl- 



■ exp 



Notice that this form for G{x) is suggested by a similar expression for the bare propagator, 

^2 



S{x) 



2nP 



exp 



x^ 
4P 



S [Xu 



(32) 



(33) 



with 



6 



... 



(see Eq. (|16D; in the chiral limit, as in the present problem, the bare mass m — Q). In order to perform the integrations 
over the perpendicular components of xi and yi in Eq. it is convenient to make use of the photon propagator in 
the momentum representation, 



(35) 



After substituting this representation along with those in Eqs. (32) and ( ^3| ) into the SD equation ( |3l| ) and performing 
the straightforward, though tedious, integrations over xj^ , y^, x^ and y^we arrive at 



g (x||) = s (x\\) + iira 



1 j2 II .2 11 

-jd x'{d y'{ exp 



+ (47ra)" 



72„ll J2^1I j2 



d'xld'yld'xld'yl exp - 



- i<l\\{x\ - 2/1 )j s(xll - x\)-i'^^g{x\ - y\)l\\g{y\)V^,,, {q\\,qi 
- ^ + ^I'ki X qi] ^qlixl - yl) - zg|(4 - y^)) 



(36) 



where [qi x q^] = €^ijqiiq2j. Since no dependence on a;^ has left, we conclude that the form of G{x) in Eq. ( p^ ) is 
indeed consistent with the structure of the SD equation. 

Regarding this equation, it is necessary to emphasize that the "perpendicular" components of the 7-matrices are 
absent in it. Indeed, because of the identity 0^~-'7'^0'-~-' — 0, all those components are killed by the projection 
operators coming from the fermion propagators. 

Substituting now the photon propagator in the Feynman-like (non-covariant) gauge ( p^ ) into the SD equation, we 

see that only the first term in Eq. (|26|), proportional to gfiu, leads to a nonvanishing contribution. In other words, the 

— . II . . . 

photon propagator is effectively proportional to gin, (justifying the name of the gauge). 

By switching to the momentum space, we obtain 



(P|l) = s ^ (Pll) - 47ra 
d^qid'^q2 



■ exp 



-(47ra)' 



(27r)8 



■ exp 



d-^q 



{qi-if 



7||'5(p" - qbil'D^ii' {q\\^q^) 



+ il'^[qi X q^ 



X7|'5(P" - 9|)7"5(P" - -zl - 91)711. 9(p" - 91)7^^- {'^Iqi) 'Dap [q'lqi) 
The general solution to this equation is given by the ansatz, 

App\\ + Bp 



g{p\\) 



Alpj-Bf 



(37) 



(38) 



where Ap = A(p|) and Bp = -B(p|). Making use of this as well as of the explicit form of the photon propagator, the 
previous equation splits into the system of two coupled equations. 



d'^qld'^q^Bp-g^Ap. 



qi-q2^p~q2 



{Ai_^,{pi\ qlY - s^,,) {AU,-qM - 9I' - qlr s^,.-,.) {Al-qAp^\ qlr Bi_ 

d^q^d^qi exp (-(q^0V2 - iqil?/^ + lP[qi X q^]) 



ql 



(g|)2n ((q^)2, (gf )2)] [q2 + (^Jl)2n ((5^X)2, (g^). 



(39) 



Bn 



d'^q\\Bp- 



'27r3 



(2713) 



(Pll-9|l) -5; 



2 

p-g 



d2g^exp(-(giOV2) 



g2 + q^ 



^(ql,qi) 



(I^q\d^ql{p\ + q\ ■ ql)Ap_q,Bp^ 



9? 



(g|)2n ((q^)2, (g^)2)] [q| + (5|)2n ((g,^)2, (q|) 



(40) 
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This concludes our derivation of the two-loop SD equations for Ap and Bp which define the fermion propagator 
through Eq. 

Let us now show that the improved rainbow (one- loop) approximation is reliable in gauge (p6|). In that approxi- 
mation, the SD equations (^9|) and ( po| ) are 

Ap = 1, (41) 



a 



d^q\\Bp^g f d^q± exp (-(qj./)V2) 



Bp = -i-^ / ^-^^^ / " ^ , ' . (42) 

Equation ( ^^ was solved both numerically and analytically (see the next section and Appendix In particular, it 
was shown that the approximation with -B(p|) = Wdyn for < 2|ei?| and i3(p|) rapidly decreasing for p| > 2|ei3| is 
a very good one (see Fig. |^ in the next Section). Moreover, as it is shown in Sec.^ and Appendix]^, in the improved 
rainbow approximation, it is sufficient to use the constant photon mass approximation for the polarization function: 
gf[n(g2^,gj|) ~ EE -2a|eS|/7r [see Eq. (|l|)]. Then, equation (|o|) at p\\ = reduces to 



a 



1 ~ — In ( , " ) In I — + Const — ^ In^ In — , (43) 

where the two terms in the right hand side come out as the estimates of the following two integrals, 

oo 


f dxidyidx2dy2 exp (-xi(TOdy„/)2/2 - X2{mdyJY /i) Jo {{mdyj)^ \/xiX2) 



87r2 J (yi + l)(j/2 + 1) {xi +yi + {M^/mdyn?) {x2 + y2 + [M^/mdynY) 



„ / yi + y2 + 1 



^{yi~y2Y + 2{yi+y2) + l 



1 . (45) 



By assuming that the first term on the right hand side in Eq. ( p3| ) is of order one, we see that the second term, 
corresponding to the two-loop correction, is indeed suppressed: it is of order alogl/a. Note that we should have 
expected this result taking into account the calculations in Appendix^. As is shown there, the higher order corrections 
to the vertex are suppressed in the Feynman-like gauge, and the two-loop contribution in the SD equation occurs as 
a result of the one- loop insertion in the vertex (see Fig. |l|) . 

If we repeat the same analysis in the case of the covariant gauge (|l^), we end up with the following estimate, 

where the main contribution to both terms, as is easy to check, comes from the '7||''z||'/9^'Z^ component in the photon 
propagator. In contrast to what we had in the Feynman-like gauge, after assuming that the first term in Eq. (^) is 
of order of one, we see that the second term would be also of order one. In other words, there is no suppression of 
higher order terms in the covariant gauge. Similarly, there is no suppression in all other gauges in which the g^g^/q^g^ 
component is non-zero. Therefore, in agreement with the general arguments of the previous section, the Feynman-like 
gauge ( |26| ) is special: in this gauge there exists a consistent truncation of the SD equations. 

Since the dynamical mass is a gauge invariant quantity, in other gauges one needs to sum up an infinite set of 
diagrams to recover the same result for it. Obviously it is rather difficult to classify all the relevant diagrams in those 



gauges, and so the existence of the Feynman-like gauge (26) is the key point that allows to treat the present problem 
reliably. 



V. THE ANALYSIS OF THE SD EQUATIONS 



As was shown in the previous section, in the gauge (p6|), the SD equations (||), (|^), (||) an d (0 ) with the bare vertex 
( p^ ) are reliable. They form a closed system of integral equations. In Euclidean space Eq. (^2^ has the form 
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a 

2^ 



dxexp{—xP/2) 



(47) 



where the polarization function IlE{q'j_,q'^^) is defined from Eq. (^) with the bare vertex ([T^). 

Equation (^) was solved by using both analytical and numerical methods. The analytical approach is considered 
in Appendix ^ Here we will describe the numerical solution. 

The results of the numerical analysis of Eq. (^^ are shown in Figs. |2[ ^ and ^. The crucial points in the analysis 
are the following: 



) is in principle a complicated functional of the fermion mass function B{p 



a) The polarization function IlE{qj^, q 
However, as it is shown in Appendix^, the leading singularity, l/Q:ln(Q!), in ln(m2^^) in Eq. (|^) is induced in the 
kinematic region with m^^,^ <C <C |ei?| and to^j,„ ^ ^ q± W^l- that region, the fermions can be 
treated as massless, and therefore the polarization function is H£;(g^,(j|) ~ 2a\eB\/TTq'^^ — M^/q'^^ [see Eqs. ( p^ and 
(pl|)]. In other words, in this approximation, the photon propagator is a propagator of a free massive boson with 
M2 = 2a\eB\/7r. 



b) In the numerical solution, the following ansatz for hi('mdyn) was used 

^dyn 



In 



0.2 



ai Na 

= In a„ + - In — - ^ 



(48) 



For small a (0.001 < a< 0.1) and different N {1 < N < 7) the best fit was found with aq = ai 
and 05 ~ 0.58 ± 0.02 (see Fig. This fit corresponds to the expression 



a2 = 03 



04 



aln(Ci/iVa) 



(49) 



— C \J\eB\F{a) exp 

where F{a) ~ {Naf^, Ci ~ 1.82 ± 0.06 and C = V2. 

c) The numerical solution shows that the function i?(p|) is essentially constant for p| ^ |e-B|, -B(p|) = mdyn, 
and rapidly decreases for ^ \^B\, see Fig. |3[ Therefore this approximation is self-consistent: the Ward-Takahashi 
identity for the vertex is satisfied in the relevant kinematic region of momenta (see Appendix |c|) , and the pole of the 
fermion propagator appears at p| = rn^dyn- 

It is instructive to clarify the origin of the difference between the dynamical mass mdyn (0), calculated in the 
rainbow approximation, and the expression for mdyn (49). The integral equation ( |47| ) for -B(p|) can be rewritten in 
the form. 



(fqB{q^ 



B^q^)J x + iq 



dx exp(— x;2/2) 

)2 + M2 ' 



(50) 



with Af^= and ^ 2a\eB\/ TT in the rainbow approximation and in the improved rainbow approximation, in 
gauge (|2q), respectively (henceforth we will omit the symbol ||; the two-dimensional vector q = (54, 53), 94 = —iqo)- 
The numerical analysis of this integral equation shows that the so called linearized approximation, with B^{q^) — s- 



'dyn 



B^ (0) in the denominator of Eq. (p(j) , is an excellent approximation. Then we get 



a f d^qB{q^ 



dxexp{—xP /2) 

.,^„, x + (q-p)2 + M2- 




(51) 



This equation is equivalent to a two-dimensional Schrodinger-like differential equation. Indeed, introducing the func- 
tion 



*(r) 



d^q B{q^) 
(27r)2 ^ 



(52) 



dyn 



we get the equation 



^dyn 



V{v))^{v) = 0, 



(53) 
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where A = ff^ /dx\ + /dx\ and the potential V^(r) is 



dx eiq){—x/2) 



X + PM^ 



(54) 



At 

Ei{x) = 



(the rainbow approximation), the last integral is expressed through the integral exponential function 
- dtexp{-t)/t, 



Voir) = Vir) 



a 



exp 



Using the asymptotic relations for Ei{x) we get 

^o(r) : 
Vo{v) : 



2a 1 
— lo 



r 

2P 



r 
'2P 



2P 



(55) 

(56a) 
(56b) 



Equation (56a) implies that the potential is long-range in the rainbow approximation. 
Using now the asymptotic relations for ^^0(2:) [p2|, we get the following asymptotes for V{r) at ^ 0, 



V{r) 



2 M 



1/2 



7 



^ I- ^log— ' 



7f 7'3/2 

^2 



r > 



> I, 



(57a) 
(57b) 



Therefore, while the short-distance behavior of the potential is independent of M-^ (compare Eqs. (56b) and (57b)), 
its long-distance behavior in the rainbow approximation and in the improved rainbow approximation is essentially 



different: while the former potential is long-range [see Eq. (56a)], the latter is short-range [see Eq. (57a)]. This point 
yields the physical origin of the difference of expressions (|l|)lHd (H) for mdyn in those two approximations. 

It is also instructive to discuss this point using results proved in the literature for the two-dimensional [d — 2) 
Schrodinger equation [ ^ . The form of equation ( [53| ) implies that —rn^ plays the role of energy E in the Schrodinger 



equation. The results of Ref. 



i.e., there is at least one solution with rn^y^ 



ensure that a) there is at least one bound state for any attractive potential for rf = 2, 
-E > 0;h) for short-range potentials of the form ^(r) = aV{r), where 



-ni 



dyn 



of the lowest bound state is —E{a) 



„(a) exp(-l/aa). 



V{r) is independent of a, the energy E 

with a > 0, as a — !■ 0. If the mass Mj were independent of a, the potential (57a) would satisfy the constraints in 
item b) above. However, because = 2Na\eB\/7r ^ as a ^ 0, we get an additional logarithmic factor in the 



power of the exponent in equation (49). On the other hand, since the potential V(r) is long-range in the rainbow 
approximation [see Eq. ( p36a| )], it clearly does not satisfy the constraints in item b) above, and rri^y^ ~ exp(— 1/a-y/a) 
in that approximation]^ 



VI. CONCLUSION 



The magnetic catalysis of chiral symmetry breaking in QED is a phenomenon with rather rich and sophisticated 
dynamics. It yields a (first, to the best of our knowledge) example in which there exists a consistent truncation of 



^As follows from the analysis in Appendix^ "^dyn would be indeed mdy„ ~ -y^jeBj exp( — 1/aa), with a ~ log |e_B|/A'/^, in 
the case of independent of a. Notice that the transition from — 2A^Q|e_B|/7r (the improved rainbow approximation) 
to = (the rainbow approximation) corresponds to changing a ~ log\eB\/M^ ~ logl/Na to a ~ log |e_B|/m^j^„. Then, 
expression (0) for mdyn transforms into expression (P). 
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the Schwinger-Dyson equations in the problem of dynamical symmetry breaking in a (3+l)-dimensional gauge theory 
without fundamental scalar fields. 

It is instructive to compare this problem with dynamical symmetry breaking in (3+l)-dimensional gauge theories 
without external fields (for a review, see Refs. [^,0). In the case of non-Abelian gauge theories, such as QCD, 
dynamical chiral symmetry breaking is generated in the infrared region, where the effective coupling constant is 
strong. This prevents to elaborate a consistent truncation of the Schwinger-Dyson equations in those theories. In 
Abelian gauge theories, on the other hand, a solution with a non-zero fermion mass exists in ladder approximation 
for any value of the coupling constant |p8| , if there is no ultraviolet cutoff in the Schwinger-Dyson equations (at finite 
cutoff, in the ladder approximation, dynamical chiral symmetry breaking takes place only if the coupling constant is 
large enough (2^]). This fact implies that the ultraviolet region is responsible for chiral symmetry breaking in that 
case. Since the running coupling is strong in the ultraviolet region in Abelian gauge theories, this again prevents to 
elaborate a consistent truncation of the Schwinger-Dyson equations in those theories either. In contrast, the dynamics 
of chiral symmetry breaking in QED in a magnetic field is long-range, and the QED coupling constant is small in 
infrared. As a result, a consistent truncation of the SD equations exists in the present problem in the gauge (|2^). 

The crucial point in the present analysis is the dimensional reduction in the dynamics of the fermion pairing. 
However, there is an essential difference between QED and the Nambu-Jona-Lasinio model in a magnetic field: in 
QED, there is an additional neutral field, the photon field. As a result, the dynamics in QED in a magnetic field 
is much more sophisticated than that in the NJL model. In particular, the photon propagator includes the four- 
dimensional = 9^ ^ (zl E^nd is not reduced to the two-dimensional form. However, the tensor and the spinor 
structure in this model in the infrared region is exactly the same as in the Schwinger model. This point is crucial for 
making this problem essentially soluble. 

What is the chiral symmetry, SUl{N) x SUr{N) x C/y(l) or Ul{N) x Ur{N), in this problem? It is known that in 
massless QED, without external fields, the chiral symmetry is Ul{N) x Uii{N): though the singlet axial current is 
not conserved, the corresponding charge is conserved. The latter is connected with the absence of instanton-like 
configurations in QED. In the presence of an external magnetic field, the situation however might be different. Indeed, 
the dynamics in QED in a magnetic field is intimately connected with the Schwinger model where the 1/4(1) symmetry 
is explicitly broken. We believe that this issue deserves further study. 

An important ingredient of the dynamics in QED in a magnetic field is the pseudo-Higgs effect. It is not a genuine 
Higgs effect since there is no complete screening of the electric charge: Eq. ( |l^ ) implies that in the deep infrared region 
with |q^| ^ rn'^y^, there are ordinary Coulomb-like forces. Still, the pseudo-Higgs effect is manifested in creating a 

resonance in the photon channel, with = 2a\eB\/Tr ^ ^"dyn- This resonance provides the dominating forces 
leading to chiral symmetry breaking. 

Thus, in this problem, the region primarily responsible for chiral symmetry breaking is the region of intermediate 
momenta: m^^^^ ^ jq^j ^ |ei3| and m^^^ ^ A/^ < \qW ^ \eB\. This point is noticeable as an example for a 
possibility discussed for the QCD dynamics [ p6| : the dynamics of spontaneous chiral symmetry breaking might be 
provided by forces essentially independent of the dynamics of confinement (the infrared dynamics) . 

Another noticeable point is the existence of a special gauge in which the description of the nonperturbative dynamics 
is essentially simplified. Recall that the conventional viewpoint now is that there is a particular gauge in QCD (the 
maximal Abelian gauge) which is the most appropriate one for the description of confinement and chiral symmetry 
breaking |29|. 

At last, there are arguments in support of a dimensional reduction in the dynamics of chiral symmetry breaking in 
QCD p3[ . The present model yields an example of such a mechanism. 

We hope that the dynamics of the magnetic catalysis of chiral symmetry breaking in QED will provide insight into 
the non-perturbative dynamics of more complicated theories, such as quantum chromodynamics. 
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APPENDIX A: LOOP CORRECTIONS TO THE VERTEX 



In this Appendix we will calculate the one and two- loop corrections to the vertex in the gauge ( p6| ) (see Fig.^). In 
agreement with the general arguments in Sec. Ill, it will be shown that they are small: of order a and in one-loop 
and two- loop approximations, respectively. On the other hand, we will show that in covariant gauges (p^, there is a 
large, correction to the vertex already in the one-loop approximation. 

As was pointed out in Sec. Ill, the loop expansion we use is based on the effective action of Cornwall, Jackiw, and 
Tomboulis. In this expansion, the full fermion and photon propagators are used and only two-particle irreducible 
graphs have to be taken into account ||2^,^. Our aim is to show that in the gauge (^6|), the improved rainbow 
approximation is reliable, i.e., the loop corrections to the vertex, with the photon and fermion propagators calculated 
in that approximation, are small. In this connection, notice that in the improved rainbow approximation the function 
in the fermion propagator equals one [see Eq. (|4^)] and -B(p|) — mdyn for all p| < 2\eB\ (see Fig. ||). Therefore 
one can use the bare propagator ([l6| ) with m = m^j,„ in the calculations of the loop corrections to the vertex: indeed, 
as was indicated in Sec. Ill, the perturbative expansion in a potentially might be destroyed by mass singularities 
coming from infrared region, where i3(p|) — mdyn- 

After getting rid of the Schwinger phase factors according to the prescription in Eq. ([l^), we get the following 
expression for the one-loop correction to the vertex function. 



r<'^^^{x -z,y-z) = (ze)2 exp {i^x - yf Ax{z ~ y)) 7^5(^)(x - z)7^^(^)(-z " y)l"T^\v{x - y) 
(see Fig. |^). By performing the Fourier transform in both x ~ z and y — z, we arrive at 

f(i)^(p, k) = (lef 



(Al) 



(27r)i 



d^x 



(2n 



(A2) 



where we first integrated over y and fc2, and then shifted fci by p + q. 

Substituting the fermion propagator as in Eq. (^), but with m — mdyn, into the last expression, we obtain 



f(i)^(p,A:) 



a 
2^ 



d^qi^d^qxexp[-{q^+p^yiy2~{q± + k^fP/2 + lP[q^ X {p^ - k^)] + lP[p^ X k^]] 



dyn 



(911 



^ii)' 



dyn 



X 7^(g|| + mdyn)0^ ^7r('7|| + k\\ + 

)0( '^Y'Dxuiq\\,q±), 



(A3) 



where [q-^ x p-^] = e^ijqiPj (with i,j = 1,2). As is clear from this expression, the dependence on the perpendicular 
momenta is weak in infrared, < |p_l|, ^ l/l, and exponentially suppressed in ultraviolet. For our purposes, it is 
sufficient to use \p±\, \k± \ — 0. Regarding the parallel components, we keep them non-zero so far. 

After substituting the photon propagator (p6|) into Eq. (|A3|), we notice that the result contains two kinds of terms, 
namely, terms proportional to 0^+' = (l-|-i7^y^sign(ei?))/2 and terms proportional to 0^~) = (1 — ?7^7^sign)(ei3)/2. 
As we discuss in Sec. [V, the former are completely irrelevant for the SD equation and we drop them. The latter are 
relevant and has to be carefully analyzed. As is easy to check, they come exclusively from the gf^, term in the photon 
propagator, and the explicit expression (in Euclidean space) reads 





7y'(9||+P||) + (a||+^||)7f 


O(-) 


in +P\\y+T^lyn 









(A4) 



In a standard way, we introduce the Feynman parameters and perform the integration over . The result reads. 



(P||, fc||) ~ -mdynO^ ^ 



oo 



X j dze ""'^ j '^^ j 




dy- 



(1 - 2/)fc||7|^ - yi\\h - ^P\n\\ + (1 - ^hiiPw 



^ipj + ^lyn) + yikn + mj ) - (xpji + yk^^)^ + {l-x~y)iz + M^) 



(A5) 
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Note that this is zero for = p|| = 0, and that the integral in the right hand side is finite for any finite rn?^y^. In 



the most interesting region of momenta, m? ^ Pf'^jj ^ l/^^i this integral contains the logarithmic contribution, 



While estimating the effect of the one- loop correction in Eq. ( A5) to the SD equation in Sec. IV, we could check that 
the most important contribution comes from the region of momenta (up to the exchange of p[| and k\\) m? <C p| <C 1//^ 
and fc|| ~ 0. In this particular case, from Eq. (A5) we obtain 

OO 2 

dze~^'' 



P\\^\\ 



■In 



—nidynO'^ '> 



'^i^m^in 



PV\\ 



In 



In 



z + M^+pI m2^„(z + M2)^ 
1 \ 



{M-,IY p\ {mdynlf {M^lf + 



(A6) 



This, in its turn, results in a suppressed two-loop correction to the SD equation. 

Now, let us consider the correction to vertex function at two-loop order. There are two diagrams in this order (see 
Fig. H). Here we shall explicitly describe the correction connected with the diagram with two crossed photon lines 
(Fig.^b); the analysis of another correction, connected with the diagram with two parallel photon lines (Fig. ^) can 
be done similarly. 

In coordinate space, the expression corresponding to the diagram with two crossed photon lines reads 



f(2)'^(x,y) = [lef j d^xid^vi exp [ix^A^ixi) + ly^A^iy) + iv^A^{x)\ -f'^ ~S^'^\x - xi)-i'' ~S^'^\xi) 

In the momentum space, this becomes 

f(^)-(p,fc) = {^er [ d'xd'yd'x^d'y^'^'^''^'^''^'^''^'^''^'''''^'''' 
X e 



(A7) 



(27r)24 

ixp~iyk+ixA(xi)+iyiA{y) — ixA(y)-iki(x-xi) — ik2Xi+ik3yi—iki(yi~y)+iqi (xi—y)+iq2{x~yi) 



X 7"5(^)(fci)7'^5(^)(fc2)7^^*^'(fc3)7''^'''nfc4)7^2?.A(gi)2?.p(g2). 



(A8) 



Performing the straightforward integrations over xi and j/i, and then over ki (by making use of the explicit form of the 
fermion propagator) , at the end we arrive at the following expression for the two-loop vertex correction at fc = p = 
(the case fc,p 7^ is not expected to give a very different estimate, but it is much harder to work with), 

f(2)M(0 0) ^ - /■ <^^<^^d^'^ld^1^d^l2 [-Wil? - Wil? + 'ill^bi X qi] 2l\t ■ q^] 



((qf )2 - m2^^„) ((g|)2 - m2^^„) (q^ + q^y - m2^„ 
X7'^(9| + mdyn)0^'h''{4 + 4+ mdyn)0^-h'"{ql + 4+ mdyn)0^~^ 



-xY{ql+mdyn)0^-~h^'D.x{qlqt)'Dap{q'^,q2) 



(A9) 



Switching to the Euclidean space and substituting the photon propagator in the Feynman like gauge (^61), we find 
that the only nonzero contribution is proportional to 0^~\ 

f(2)A'(0,0) = ^^^^7''O(-) 

d^q.d^q^d^q^d^qi exp [-|(g^0' - + ^^I'kt X qi] - 2Pqt ■ q^] 



qf + ^dyn [qi + ^dy 



(qi + 92) 



'■dyn 



{qf + {qir + M^) {ql + {qir + M^) 



< Const7'fO(-)a2 In^ 



1 + {M^iy 



(AlO) 
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If we assume that N is not too large, so that {M^l)^ ~ 2a/7r ^ 1, then the above estimate for the two-loop correction 
becomes 

f (2)M(o, 0) < Const 7^0(-)a2 W (^) . (All) 

Thus, we conclude that, as expected, the Feynman like non-local gauge leads to suppressed higher order corrections. 
The latter, in its turn, means that the solution to the SD equation with the bare vertex and the photon propagator 
in Eq. (|2^) is reliable, and that the result presumably approaches the exact one when the fine structure constant is 
very small. 

At this point it is instructive to explicitly demonstrate that the observed suppression in higher orders is the exclusive 
property of the special gauge in Eq. (26). 

To prove this, let us consider th e on e-loop vertex correction in the case of the covariant gauge (0). After substituting 



the photon propagator into Eq. (A3), we arrive at the following estimate 



f<.V(0,0).ii.fO<-.,.,'(^)-ii.fO<-.,n 



M'^ \, / 1 



In 



where, as is easy to check, the leading term with double logarithm comes from the gjigll/fj'^fj'^ component in the photon 
propagator in Eq. (|l^). Since in the rainbow approximation the solution to the SD eq uation in the covariant gauge 
yields \ia^{mdyniy ~ 1/a [see Eq. (0)], we conclude that the one-loop correction in Eq. (A12) is of the same order as 



the bare vertex, i.e., f ^^^''(0,0) — 0(1). As a result, such an approximation is not self-consistent. 

Notice that if were equal to zero (i.e., Il{q'j_, q'^^) = in Eq. ( |l7| ) and the gjlgl/g'^Q^ term would be absent there), 
the first two terms in Eq. ( A12| ) would cancel, and the one- loop correction f(i)'^(0,0) would be small (f(i)^(0,0) 



0{^/a)). As a result, expression (^ for rudyn would be correct. Therefore the origin of the deviation of expressions 
(|l|) and (||) for nidyn is the generation of a non-zero M^. 

We emphasize that the qjiql/q^q^]^ component in the photon propagator, responsible for breaking the consistency 
of the rainbow approximation to the SD equation, is absent in the gauge Eq. (|6|). As a result, in this gauge, there 
exists a consistent truncation of the SD equation for small values of the coupling constant. Since the dynamical mass 
is a gauge invariant quantity, in other gauges one needs to sum up an infinite set of diagrams to recover the same 
result. Obviously it is rather difficult to classify all the relevant diagrams in those gauges, and so the existence of the 
special gauge is the key point that allows to solve the problem of the magnetic catalysis in QED reliably. 



APPENDIX B: ANALYTICAL ANALYSIS OF THE SD EQUATION 

In this appendix we will describe the analytical solution of the SD equation (^). 

First of all, let us show that the leading singularity, l/aloga, in log(TO^j^„) in Eq. ( |49| ) is induced in the kinematic 
region with m\^^ <^ <C \eB\ and rn^y^ <^ ^ <C \eB\ (in that region, fermions can be treated as massless). 
As was shown in Sec. the approximation with i3(p|) = mdyn for p| < 2|e_B| and -B(p|) rapidly decreasing for 
> 2\eB\ is reliable in this problem (see Fig. |3|). Then, taking p| = in Eq. (^7|), we arrive at the equation 

^2 J dxeiip(-xP/2) 



a 



a 



''dyn 
'II 



27r2 J qf^+mj J x + qf^+qf:nE{x,qf:) 







, , 2|eB| 

d^q« f dx 



2n^ J qfi+mly^ J x + q^^+qjnE{x,qp' 
Matching now the asymptotes (18),(p^) at (/jj = 6mdyn in Euclidean space, we get 



(Bl) 



14 



2\eB\ 



a 



2\eB\ 



dx 



6m^ 



dy 



i iy + mXj{x + y{l + ^) 



dy 



dyn 

It is clear that, because of in (y + n^dyn)' ^^'^ ^^'^^ term in the square bracket on the right hand side of this 

equation is of order 0(1) and can be neglected: it cannot give a contribution of order 1/aloga to Then we 

arrive at the estimate, 



(B2) 



1^^ 



2\eB\ 2\eB\ 

dy 



2n J y + m\y^ 

6mj 



dx 



x + y + M2e~=^'V2 ' 



(B3) 



The double logarithmic contribution comes from the region 2\eB\ ;:|> y = ^ ''^^dym 2|eB| ^ a; = > y + > 
M^, = 2a /nP. Therefore one can write 



2\eB\ 2\eB\ 




2\eB 


a r dy r 


dx 


a r 


27T J y J 


X 


27r J 






6m ^ 

dyn 



log- 



(B4) 



To calculate the last integral with double logarithmic accuracy, we write 

2\eB\ 



2-K 



a 
2^ 



log 



2eB 



dy 

y 



dy, 2\eB\ 
— log' 



6mi 



Mi 



y y 



, 2eB , M2 1 2eB 
T72 log + o T72 



a , 2eB , 



2\eB\M^ 



This equation implies that 



^dyn 



Na 

TT 



1/4 



exp 



(B5) 



(B6) 



Comparing this expression with Eq. (p9|), one can see that this estimate is quite reasonable. The origin of that is a 
rather simple form of the fcrmion mass function — rudyn for p| ^ 2\eB\ and rapidly decreases for 

pjj >2|ei?|. 

Therefore the dominant contribution to the SD equation (p7| ) comes from the region with nn\^^ <C |q^| <C \eB\ and 
Tn^dyn -^^7 ^ <fi ^ k-^li whcrc fermions can be treated as massless. This in turn justifies the approximation with 
the polarization function He ~ 2a|ei3|/7rg^. 

Now we proceed at solving analytically SD equation ( p7| ) for a mass function in the improved ladder approximation. 
It is 



27r2 y k^ + B^{k)J z + {k-pY + T\{z)'' 



(B7) 



where = 2a\eB\/TT and we shifted the momentum of integrati on (also we omit the symbol || in the rest of this 
appendix). After the integration over the angular coordinate, Eq. (B7) becomes 



a f dk^B{k^'' 



K{p\k^) 



k^+B^{k^) 



(B8) 
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with the kernel 



K{p\e) 



dz exp(— 



(B9) 



To study Eq. (B8) analyticaUy, we break up the momentum integration into two regions and expand the kernel 
appropriately for each region (compare with Refs. pj 



dz exp(— 0^/2) 



k^ + B^{k^)J p2 + ^2g-2/V2 + z 



dPB{k^) f dzeM~zP/2) 



k2 + B^k^)J k^ + IVne'-^y^ + z 



(BIO) 



Introducing dimensionless variables x = y — k'^P /2 and also the dimensionless mass function B{p'^)/ ^/2\eB\ 

B(x), we rewrite the last equation in the form 



B{x) 



2ti 



dyB{y) f dyB{y)g{y) 



y + B^{y) J y + B^y) 



where 



(Bll) 



dze 



z + X ^ — e 



(B12) 



The solutions of the integral equation ( Bll ) satisfy the second-order differential equation 

B 



B"-3—B'~—q'- 



= 0, 



g' 27r-^ x^B'^ix) 

where the prime denotes derivative with respect to x. The boundary conditions are 



0, 



(B13) 



(B14) 



B 



gB' 



= 0. 



The function g{x) has asymptotic behavior 



(B15) 



9KX) 



log- 



X < 1, 



a; > 1. 



(B16) 



We consider now the linearized version of Eq. (B13) when the term B'^{x) in denominator is replaced by a constant 
_B2(0) = a? {B{p = 0) = radyn)'- the numerical analysis shows that it is an excellent approximation. The two 
independent solutions of that equation near the point x = oo behave as B{x) const and B{x) l/x, and the 
UVBC selects the last one. 

In the region x <C 1, the equation takes the form 



B" 



Kb' " 



B 



2^ (x+f)(x + a2) 



0. 



(B17) 
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Introducing the variable x + = —ziajir — Eq. (B17) can be rewritten in the form of an equation for the 
hypergeometric function, 



z(l-z) 



£-B dB a 



dz^ dz 27r 



B = 0. 



The general solution to Eq. (B18) has the form 

B{z) = Ciui + C2U2, 

where 

Ml — zF{l + — iv] 2; z), 

U2 = {-z)-'''F{w, l + iv]l + 2iv; ^) + {~zy''F{-iv, l-iv;l- 2iv; i), 
V — yJa/2'K. From the infrared boundary condition ( [B14|) , one gets 

Ci u' 

^ x=0 



(B18) 

(B19) 

(B20) 
(B21) 

(B22) 



Equating now the logarithmic derivatives of solution (B19) (at a; <C 1) and 1/x (at a; 3> 1) at the point x — 1, we 
arrive at the equation determining the quantity a (a) {i.e., the dynamical mass nidyn)- 



ip = A1B2 ~ A2B1 = 0, 



where 



Since the variable z is 



A, = (u^ + u,) 



B, = 



x=0 



X + a'^ 



X + a'^ 



-—a , z 
a 



a 



(B23) 



(B24) 



(B25) 



(we suppose ^ ct/ir), in what follows we need asymptotic behavior of Ui{z), u'^{z) at small and large negative values 
of its argument z. Using corresponding formulas from |]30|] , we find for small values of z {\z\ <^ 1) 



7.1 ~ z(l + i±^z) + 0(z2) 



U2 ~ 2Re 



u[~l + {l + iy^)z + 0{z^), 



u'2 ^ 2Re 



where 



^^^^„^i\ni-z) + l-ho) + 0{z\nz) } , 
/lo = 1 - 27 - ^{w) - tp(l + if). 



At large \z\ >• 1 we have 



Ui 



1 / tanh(7ri/) 



sin [z/ln(-z) + $(i/)] + 0(z"i), 



U2 ~ 2cos(iyln(-z)) + 0(z"^), 



1 / tanh(7ri') 



2 V TTiy 
4 ^ -^sin(i^ln(-z)) + 0(z"^), 



cos(i'ln(-2;) + $(i^)) +0(z 



(B26) 

(B27) 
(B28) 
(B29) 

(B30) 

(B31) 
(B32) 

(B33) 

(B34) 
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where 



^(v) = are 



T{l + 2iv) 



_i2(22"-l)C(2n+l) 



,2n+l 



r^{l + iiy)J ^' 2n + l 

~ 2C{3)iy^ - 6C(5)i^^ + . . . . 
By making use of these asymptotes, we obtain the following expressions for Ai, Bi 



tanh(7ri/) 



cos(iy m — ) H — sm I m — 
a \ a 



A2 
B2 



,dU2 „ r,-/i''''\ 

— ; h M2) U=i — 2 cos \vm—] +2 sm U^ln — 1 , 

dx \ aJ \ aJ 

dui 



dx 
du2 
dx 



\x=0 



\x=0 



TT 

a 

-In- 



And, finally, the solution to Eq. (B23) reads 



2 _ "^dyn 



Na 



2\eB\ 



■ exp 



1 / , TT 

- — cot 1/ In — — 
V \ Na 



2/3 



exp 



27r 



a ln(7r/A'^Q!) 



as a ^ 0. 



(B35) 
(B36) 

(B37) 

(B38) 
(B39) 
(B40) 

(B41) 



The obtained analytical expression (B41) for the dynamical mass is close both to estimate (p6|) and to the numerical 



(anal) />^(numcr) 



1.7. 



solution (^9|). The ratio of the values of Ci in the analytical solution and in the numerical one is jC- 
This, rather mild, discrepancy reflects the approximations made in the kernel of the integral equation (|B8) in reducing 
it to the differential equation (B13). 



APPENDIX C: WARD-TAKAHASHI IDENTITY 



In this appendix, we discuss the simplest Ward-Takahashi (WT) identity relating the fermion propagator and the 
vertex in the problem of magnetic catalysis (for another approach see Ref. p^). 
Let us start from the WT identity in coordinate space: 



d^^V^{x', y'- z) = 8{x' - z)G-\z, y') - 6{z - y')G-\x' , z) 



(CI) 



After multiplying this expression with G{x,x') on the left and G(y', y) on the right and integrating over x' and we 
arrive at another representation. 



d^x'd^y'G[x, x')dlT^{x\ y'- z)G{y\ y) - 5{z - y)G{x, z) - 5{x - z)G{z, y). 



(C2) 



By taking into account the universal Schwinger phase factors of the propagators and the vertex [see Eq. (|T^) in Sec. 
this identity (in the momentum representation) reads 



iG{p2) - iGipi) = / dV^e— «G (pi + A^-^\r)) (p, - pz)^? ^(pi + q,P2 + q)G {p2 + q - A^^*(r)) 



(C3) 



Now, by substituting the bare vertex along with the solution for the fermion propagator as in Eq. (|28| ) with j4(p|) = 1 
and -B(p|) = rndyn = Const (recall that the latter is a very good approximation in the most important region of 
momenta, y <C \eB\), we obtain the following relation 



(pf - mdyn)e^ 



(piif _ (J _ 



(p| - m,,„)e-(^'^') O(-) = e-(ptly/2-(ptly/2+^l^lPt -Pt] (pf _ . 



(C4) 



Once again, the approximation for the fermion propagator with A(p^) — 1 and i?(p^) — rndyn = Const is reliable 
only in the dynamical region ,, ^ \eB\. Therefore, at best, one could expect that the WT identity is satisfied in 
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this same region of momenta. After expanding b oth sides of Eq. ( |C4| ) in powers of p]_P = p\/\eB\ and keeping only 
the leading order contribution, we see that Eq. (C4) indeed turns into an identity. In the end, we remind that the 
same conclusion about the WT identity was reached in Ref. by making use of a different approach. 
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FIG. 1. The diagrammatic representation of the SD equation in coordinate space at two- loop order. A thin solid line 
corresponds to the bare fermion propagator; a bold line corresponds to the full fermion propagator, and a wavy line designs 
the full photon propagator. 
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FIG. 2. Plot of function of a for several values of Nf . 
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FIG. 3. Plot of the mass function B{p^) as a function of for Nf = 1 and two values of a. 
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FIG. 4. Plot of the fit function and corresponding datapoints for several values of Nf when the only free parameter was os- 
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